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The probability distribution of log-returns for financial time series, sampled at high frequency, is 
the basis for any further developments in quantitative finance. In this letter, we present experimental 
results based on a large set of time series on futures. We show that the t-distribution with i/ ~ 3 
gives a nice description of almost all data series considered for a time scale At below 1 hour. For 
At > 8 hours, the Gaussian regime is reached. A particular focus has been put on the DAX and Euro 
futures. This appears to be a quite general result that stays robust on a large set of futures, but not 
on any data sets. In this sense, this is not universal. A technique using factorial moments defined 
on a sequence of returns is described and similar results for time scales are obtained. Let us note 
that from a fundamental point of view, there is no clear reason why DAX and Euro futures should 
present similar behavior with respect to their return distributions. Both are complex markets where 
many internal and external factors interact at each instant to determine the transaction price. These 
factors are certainly different for an index on a change parity (Euro) and an index on stocks (DAX). 
Thus, this is striking that we can identify universal statistical features in price fluctuations of these 
markets. This is really the advantage of micro-structure analysis to prompt unified approaches of 
different kinds of markets. Finally, we examine the relation of power law distribution of returns 
with another scaling behavior of the data encoded into the Hurst exponent. We have obtained 
H = 0.54 ± 0.04 for DAX and H = 0.51 ± 0.03 for Euro futures. 

Keywords: quantitative finance, distribution of returns, t-distribution, non-extensive statistics, factorial 
moments 



I. INTRODUCTION 

Returns in financial time series are fundamental inputs to quantitative finance. To a certain extend, they provide 
some insights in the the dynamical content of the market. There are several experimental analysis showing that the 
probability distributions of returns in a large set of financial markets exhibit power law tails [iH^. 

Let us consider financial price series, labeled as S{t), from which we extract the log-returns x{t) — log[S'(t-|-A<)/5'(t)] 
over some time interval At. Any statistical analysis can then be conducted on these log-returns x{t). There are 
evidences that for mature and hi gh l iquid markets, in particular futures, both positive and negative tails conform to 
the so-called inverse cubic law . It means that if we express the the distribution of returns as a power law 

/(H) (1) 

at large values of |a;|, we can measure the exponent ^ ~ 3 as it is done in Ref. [^, [l3- This is correct for a variety of 
mature markets, in particular futures. This observation is in contrast to p redictions from the Pareto-Levy distribution, 
and former experimental results, from which we expect 1 < i/ < 2 [ll|. This discrepancy in results shows how the 
topic is critical and a clear driver for theoretical considerations. There is no universality and the use of large samples 
of high frequency data is a key point in the understanding of the underlying market dynamics. 

Very generally, all experimental investigations confirm that distribution of returns evolve from power law behavior 
at small time scales At (see Eq. ([1])) to Gaussian at large time scales. The precise value of the exponent i/ depends on 
the market under consideration with a tendency for mature and liquid markets to fall outside the Levy stable regime, 
namely v > 2. 

In this letter, we consider a few financial series on futures, which correspond to mature markets with high level 
of liquidity. First, we discuss the universality of the experimental distributions of returns on these financial series. 
Following Ref. [lOj, we extract a parameterization of the log-returns distribution which confirms the previous result, 
j/ ~ 3, for small time scales At. 

Then, we show at which value of the time scale, market dynamics enters into the Gaussian regime. We discuss the 
universality of this particular scale with respect to a sample of high liquid futures. In a second part, we address the 
determination of these scales using an analysis technique based on factorial moments [T^ . We prove the consistency 
between the appearance of intermittency and the deviation from the Gaussian regime, as derived in the first part. 

From the measurement of the distributions for returns x{t) = log[5'(t -b At)/S{t)], we can understand how these 
returns sum up to build the quantity x{t, r) = log[5'(t -I- r)/S'(t)], where t is a multiple of At. Then, in a thirst part, 
we discuss the scaling of < x{t,T)‘^ >t with r. We relate the so-called Hurst exponent to the power law tail i/. 
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II. DISTRIBUTION OF RETURNS P{x)[At] 

A. Definitions 

From standard quantitative analysis we know that the distribution of log-returns x{t), namely P{x), can be 

written quite generally as 

P{x) = ^exp{-^w{x)/2) (2) 

where w{x) is an objective function and Z a normalization factor. In particular, it can be shown easily that w{x) 
can be derived by minimizing a generating functional F[w{x)], subject to some constraints on the mean value of the 
objective function. In Eq. ([2|) we do not specify the time scale At at which returns are derived. When need, we will 
take this variable into account with the notation P{x)[At], It reads 

F = y dxP{x)[\ogP{x) +w{x)/D — X] (3) 



where A is an arbitrary constant. 

In addition, the expression given in Eq. © for the probability distribution can also be seen as the outcome of an 
equation of motion for x{t). From Eq. ([2]) and ([3]), we can express the stochastic process x{t) as a Markovian process 
of the form 



= f{x) + g{x)e{t) (4) 

where e{t) is a Gaussian process satisfying < e{t)e{t') >= DS{t — t') and < e{t) >= 0. In Eq. (|3]), functions / and 
g depends only on x{t). Adopting the Ito convention M, the distribution function P{x,t), associated with this 
equation of motion (Eq. Q), is given by the following Fokker-Planck equation 

dPix,t) ,D ,, . , ,, 5 , . , , ,, 

- ^[nx)P{x,t)] (5) 

From Eq. we can finally extract the stationary solution for P{x) in the form of Eq. ([2]) 



exp 
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( 6 ) 



Whether Eq. and (HI) can be related to real data on financial markets is not granted. Therefore, we need 

to compare predictions derived from these equations to real data. As mentioned above, we use financial time series 
on different futures, using a five minutes sampling. 

In Eq. (H]), dS|) and dH), functions / and g are not specified and any choice can be considered. Obviously, only 
specific choices will have a chance to get a reasonable agreement with real data. For example, let us consider three 
cases: 



If f{x) = —X and g{x) = 1, we obtain P{x)Z = exp(— x^/U), and thus we predict a Gaussian shape for the 
log-returns distribution. 

In the more general case where /(x) = Xg^ and g is not constant, we obtain 

= g2(l-\/D) 

and thus we predict non-Gaussian shape for the log-returns distribution. 

(iii) Let us specify the case (ii). Introducing a constant zz and defining the two functions / and g as /(x) = 
and g{x) = (1 + ^), we get 

P{x)Z = — — . 

^ ^ (1 + ^)( i ^+ l )/2 



(i) 

(ii) 



( 7 ) 
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In this scenario P{x) follows the so-called t-distribution. It depends on one parameter i/ to be fitted on real 
data, for normalized log-returns. Let us notice that Eq. © is equivalent to the q-exponential form of Ref. 



[TM3 



P{x)Z 



1 

(1+x2§5|)1/(9-1)' 



( 8 ) 



where we have conserved the notations of Ref. [l3l - [l5l| . Obviously, Eq. © and © are directly related by 
(ly -\-l)/2 = l/{q — I). In particular, j/ = 3 is equivalent to q = 1.5. 




5min Log-returns 



FIG. 1: Log-returns x (five minutes sampling) for a large set of futures. To make the comparison, we have scaled x for all 
futures to the volatility of the DAX futures (FDAX). we observe that futures on DAX, Bund, Yen, Euro, Gold present the 
same probability distribution P(x). 



B. Experimental Analysis of P{x)[At] 

In Fig. [l] we present the log-returns x (five minutes sampling) for a large set of futures, with At = 5 minutes. To 
make the comparison, we have scaled x for all futures to the volatility of the DAX futures (FDAX). All data sets 
cover the period 200I-20II. On the left hand side of Fig. [U we observe that futures on DAX, Bund, Yen, Euro, Gold 
present the same probability distribution for x. Therefore P{x) is universal for all these data series once the volatility 
is normalized to the same value. Note that there are futures on commodities which exhibit some larger tails and 
we can not claim the universality of this distribution whatever futures. See also Ref. Q. 

We can use results developed above in order to compare with experimental distributions P{x) of Fig. [T] (with fixed 
At = 5 minutes). We use predictions exposed in cases (i) and (iii), respectively the Gaussian and the q-exponential 
forms (or t-distribution). Results are shown in Fig. For data, we only display P{x) for DAX futures. 

In Fig. [5J we show that the q-exponential probability distribution of Eq. ®, with q = 1.5, gives a good description 
of the data. Similarly, it corresponds to = 3 in the form of the t-distribution of Eq. ©. Also, we observe in Fig. 
[5] that the Gaussian approximation fails to describe properly the data. On the right hand side of Fig. ^ we observe 
that when v is increased above 3 in Eq. ©, then P{x) stands in the middle of the correct probability density and 
the Gaussian approximation. When v tends to infinity, the t-distribution recovers the Gaussian limit. 
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FDAX log returns q-Gauss and t-dist 



FIG. 2: We show that the q-exponential probability distribution, with q = 1.5, gives a good description of the data (Right). It 
corresponds to = 3 in the form of the t-distribution (Left) . The Gaussian approximation fails to describe properly the data. 
When 1 / is increased, the Gaussian limit is approached. 



In summary, from the theory point of view, the probability distribution of log-returns of financial time series, 
sampled at high frequency, can be expressed quite generally as P{x) = exp{—w{x)/ D). In particular, we have 

shown that a large sample of high liquid futures (with normalized log-returns) are compatible with a distribution of 
the form 



P{x) (X 



1 

(1 + ^)( i ^+ l )/2 



~ 3. 



(9) 



Commodity futures deviate from this shape with larger tails and thus a smaller value of the exponent v. 

Once Eq. m is established, we can examine the validity of this probability distribution as a function of At, from 5 
minutes to several hours. This is done in Fig. [3] for the DAX futures (FDAX) and Fig. 0]for the Euro futures (EC). 
As mentioned above, these two data sets cover the period 2001-2011. In Fig. [3] (a)-(d), we present the normalized 
log-returns for the FDAX on different time scales (or resolutions) At. The normalization is done by changing x in 
{x — fix) / fXx where fix and ctx are the sample mean and variance. This re-definition does not change obviously the tail 
behavior of P{x). 

The basic resolution At = 5 minutes is displayed in all cases (a)-(d) as a reference. Then, we show the distributions 
using At = 10,20 minutes for Fig. [3] (a). At = 30,60 minutes for Fig. [3] (b). At = 2,4 hours for Fig. [3] (c) and 
finally At = 8, 16 hours for Fig. [3] (d). For small times scales At < 1 hour, in Fig. [3] (a)-(b), we confirm that the 
distributions of log-returns follows Eq. with z/ = 3. 

When the resolution is increased to At = 2,4 hours. Fig. [3] (c), we observe that P(a:)[At] starts approaching the 
Gaussian regime. This regime is reached for At > 8 hours as can be seen on [3] (d) . 

In Fig. |4](a)-(d), we present similar plots for the Euro futures. Same conclusions can be derived. For small times 
scales At < 1 hour, in Fig. 0] (a)-(b), we observe that the distributions of log-returns follows Eq. ([3]) with v = ?>. 
When the resolution is increased to At = 2,4 hours. Fig. 0] (c), we observe that P{x)[At] starts approaching the 
Gaussian regime, which is reached for At > 8 hours (01(d)). 

Identical conclusions are valid also for futures displayed in Fig. [T] This is a reasonable behavior which is compatible 
with the general statement done in the introduction. P(a:)[At] evolves from power law at small time scales At to 
Gaussian at large time scales, namely At > 8 hours. With the above analysis, we have put quantitative estimates on 
this statement for several high liquid futures. 
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FDAX log-returns 



FDAX log-returns 



FIG. 3: Normalized log-returns for the FDAX on different time scales At. The normalization is done by changing x in 

{x — fix)l<yx where Hx and ax are the sample mean and variance. This re-definition does not change obviously the tail behavior 
of P{x). The basic resolution At = 5 minutes is displayed in all cases (a)-(d). Then, we show the distributions using At = 10, 20 
minutes (a). At = 30,60 minutes (b), At = 2,4 hours (c) and finally At = 8,16 hours (d). For small times scales At < 1 
hour (a)-(b), we confirm that the distributions of log-returns follows Eq. ([9} with zz = 3. When the resolution is increased to 
At = 2, 4 hours (c), we observe that we start approaching the Gaussian regime. This regime is reached for At > 8 hours as can 
be seen on (d). 
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EC log-returns 



EC log-returns 





EC log-returns 



EC log-returns 



FIG. 4: Normalized log-returns for the Euro futures (EC) on different time scales At. The normalization is done by changing 
a: in (x — j Ox, where and ax are the sample mean and variance. This re-definition does not change obviously the tail 
behavior of P(x). The basic resolution At = 5 minutes is displayed in all cases (a)-(d). Then, we show the distributions using 
At = 10, 20 minutes (a). At = 30, 60 minutes (b). At = 2, 4 hours (c) and finally At = 8, 16 hours (d). For small times scales 
At < 1 hour (a)-(b), we conhrm that the distributions of log-returns follows the same distribution probability as the FDAX. 
When the resolution is increased to At = 2, 4 hours (c), we observe that we start approaching the Gaussian regime. This regime 
is reached for At > 8 hours as can be seen on (d). 
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III. FACTORIAL MOMENT ANALYSIS 
A. Definitions 



In this section, we intend to analyze the dynamics of returns using an alternative technique. The main interest is to 
show the consistency of all results within different techniques. In the following, we are interested in the multiplicity 
of positive or negative returns in a given time window At. Indeed, sequences of positive and negative returns are 
much indicative of the statistical nature of fluctuations in the price series. The idea is then to extract a quantitative 
information from these sequences. 

First, we present the situation in nuclear physics. At nuclear or sub-nuclear energies, the number of hadrons 
created during inelastic collisions varies from one event to another. The distribution of the number of produced 
hadrons, namely the multiplicity distribution, provides a basic means to characterize the events. The multiplicity 
distribution contains information about multi-particle correlations in an integrated form, providing a general and 
sensitive means to probe the dynamics of the interaction. Particle multiplicities can be studied in terms of the 
normalized factorial moments 



Fq{T) = (n(n - 1) . . . {n - q + 1))/ (n)«, q = 2, 3, . . . , 



(10) 



for a specified phase-space region of size T. The number, n, of particles is measured inside T and angled brackets 
(. . .) denote the average over all events. The factorial moments are convenient tools to characterize the multiplicity 
distributions when T becomes small. For uncorrelated particle production within F, Poisson or Gaussian statistics 
holds and Fq = 1 for all q. Correlations between particles lead to a broadening of the multiplicity distribution and to 
dynamical fluctuations. In this case, the normalized factorial moments increase with decreasing F. The idea is then 
to divide the factorial moment defined in Eq. (Hoj in more and more bins. 

We can thus compute the related moment following Eq. (HQl) as 



F, = 



1 



N. 



events 



E 

events 
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{{n)/nbinsY 



( 11 ) 



where {n) is the average number of particles in the full phase space region accepted (F), nuns denotes the number of 
bins in this region and is the multiplicity in fc-th bin. 

The behavior of factorial moments plotted as a function of the number of bins (which means decreasing bin sizes) 
provides information about the character of multiplicity fluctuations among different bins. Rising of Fi with rising 
nbins (decreasing bin size) generally signalizes deviation from purely Gaussian distribution of fluctuations. The linear 
growth of log Fi with nuns was defined as intermittency in [T^ . See also [Sm. In the following, the term is used 
for any type of growth of Fi observed. 

As a matter of fact, it has been noticed in [I^ that the use of factorial moments allows to extract the dynamical 
signal from the Poisson noise in the analysis of the multiplicity signal in high energy reactions. 

In addition, it has been shown that it is possible to define and compute a multi-fractal dimension, 1?^, for the theory 
of strong interactions [2l|, 

Fq{T) = {n{n -l)...{n-q + l))/(n)« oc (12) 

where d is the dimension of the phase space under consideration [d = 2 for the whole angular phase space, and d = 1 
if one has integrated over, say the azimuthal angle). In the constant coupling case Vq is well defined and reads 

(13) 

where 7 q = a>^s ^ cts is the strong interaction coupling constant, Ca is the gluon color factor [U The 
choice of the factorial moments as a specific tool in order to study the scaling behavior of the high energy multiplicity 
distributions is then useful to analyze the underlying dynamics of the processes. In principle, we can extend this last 
idea to other fields where factorial moments can be defined. 



B. Generating Function for Factorial Moments 

The multiplicity distribution is defined as where cr„ is the cross section of an n-particle produc- 

tion process (the so-called topological cross section) and the sum is over all possible values of n so that 

OO 

Y.Pn = l. 

n—0 



(14) 



The generating function can be defined as 



G{z) = J2Pn{l+Zr, 

n—0 



( 15 ) 



which substitutes an analytic function of 2 in place of the set of numbers Then, we obtain the factorial moment 
or order q as 



1 


d‘iG{z) 




(n)9 


dzi 


Z — 



(16) 



and the corresponding definition for cumulants 



_J_ d<i InGjz) 
{n)i dzi 



(17) 



The expression for G{z) can then be re-written as 




(18) 



\nG{z) = Y,"-{nYK, {K^ = l). (19) 

9=1 

The physical meaning of these moments has been discussed in the previous section. Another interpretation can be 
seen from the above definitions if they are presented in the form of integrals of correlation functions. Let the single 
symbol y represent all kinematic variables needed to specify the position of each particle in the phase space volume T 
|l2l |. A sequence of inclusive g-particle differential cross sections d'^a/dyi . . . dyq defines the factorial moments as 

~ 7 ^ f dyi... [ dyq ^ ^ . ( 20 ) 

{n}'^ Jr Jr dyi...dyq 

Therefore, factorial moments include all correlations within the system of particles under consideration. They represent 
integral characteristics of any correlation between the particles whereas cumulants of rank q represent genuine q- 
particle correlations not reducible to the product of lower order correlations. 



C. Experimental Analysis 



The analogy with returns of financial price series is immediate. If we divide the price series y{t) in consecutive 
time windows of lengths At, we define like this a set of events. In each window, we have a certain number of positive 
returns n+, where y(t) — y(t — I) > 0, and similarly of negative returns n_. If the sequence of returns is purely 
uncorrelated, following a Gaussian statistics at all scales, we expect Fq = 1 for all q. 

However, correlations between returns may lead to a broadening of the multiplicity distributions (n+ or n_ or 
even a combination of both) and to dynamical fluctuations. In this case, the factorial moments may increase with 
decreasing At, or increasing the number of bins that divide At, as in Eq. (fTTl) . In the following, we consider only the 
factorial moment of second order T 2 . We can write 



1 Sfc'Oi'’ (’^fc 1)} 

Prevents / 'IT’binsY 



where (n) is the average number of positive returns in the full time window (A), nbins denotes the number of bins in 
this window and n'^ is the number of positive returns in A:-th bin. 

We consider the data series on the Euro futures, sampled in 5 minutes units (same data set as in section ID). We 
present calculations for a time window 200 times units, that we divide in 1, 2, 4, 10 or 20 bins. This means that 
the time resolution extends from 1.6 hours to 16.5 hours. Note that with a time window of 200 time units, we set 
up an ensemble of more that 3400 events. The statistical precision of the following analysis is then ensured. Results 
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FIG. 5: Euro futures. Left: factorial moments are displayed for 1, 2, 4, 10 or 20 bins. The global time window of 16.5 
hours (ubins = 1). This provides a time scale ranging from 1.6 hours for 20 bins till 16.5 hours for 1 bin. We observe that for 1 
and 2 bins segmentation, is found equal to 1 and F^~^ is increasing above 1 when the number of bins gets larger than 4. 
This confirms that non-Gaussian fluctuations in the sequence of returns returns start to play a role when the scale (resolution) 
is below 4 hours. Right: similar results and comments for unlike-sign factorial moments F^ . Note that intermittency is 
enhanced for F^ compared to the like sign quantity F^~^. 



are presented in Fig. [5] (left) for positive returns. Factorial moments are displayed for 1, 2, 4, 10 or 20 bins. 
As mentioned above, this gives a time resolution ranging from 1.6 hours for 20 bins till 16.5 hours for 1 bin. The 
statistical precision is of 0.1%. We can define a systematical uncertainty by shifting the time window of 200 units 
by 50 or 100 units, which means that we define a different set of events among the price series. Variations in the 
calculations of are lower that 0.1%. Fig. [5] (left) displays the full uncertainty of these quantities. 

In Fig. O (left), we observe that for 1 and 2 bins segmentation F^~^ is found to be equal to 1. As expected, for 
the larger resolution, positive returns appear as completely uncorrelated. When the number of bins is increased, we 
observe the phenomenon described in section 1, with an enhanced sensitivity of ^ 2 *”^ to non-Gaussian fluctuations. 
This confirms that correlations between positive returns start to play a role when the resolution is below 4 hours. 
Thus, Fig. [S] (left) exhibits a clear feature of intermittency. As in nuclear reactions, an increase of the resolution to a 
certain extend leads to a broadening of the multiplicity distribution and to super-diffusive fluctuations. Let us note 
that this is a feature that can be approached in the context of non-extensive statistics [isl - fT^ . 

Similar results can be obtained for F 2 , defined for negative returns distribution. 



FF~ = 



w. 



{^fc 1)} 

{{n~) / UbinsY 



( 22 ) 



where (n) is the average number of negative returns in the full time window (A), Ubins denotes the number of bins in 
this window and is the number of negative returns in fc-th bin. For all values displayed in Fig. [S] (left) for Fj*”'', 
we derive the same result for F 2 up to 0.1%, which makes and F 2 indistinguishable. 

A direct extension of the above study can be obtained if we examine moments for like-sign and unlike-sign combi- 
nations of returns separately. The like-sign factorial moment of order 2 is defined by Eq. (EH- The unlike-sign can 
be expressed as 



F'' — 



W, 
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events 






(23) 
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Results are presented in Fig. [5] (right) for F 2 . Here again, we observe intermittency, with an increase of F 2 as a 
function of the number of bins. Also, this increase is larger than for like-sign calculations as can be observed when 
comparing Fig. [5] (left) and [5] (right). The sensitivity to non-Gaussian fluctuations in the returns sequence is thus 
enhanced with the definition of d to F 2 . This is also an effect observed in nuclear interactions . 

: 1.2 

CM 

“- 1.175 
1.15 
1.125 
1.1 
1.075 
1.05 
1.025 
1 

0 5 10 15 20 25 

nb of bins 



II II 

^ E 
0 F 


1 1 1 1 1 1 1 1 

iurofuti 

■DAX 


1 1 1 1 

ire 

C 


1 1 1 1 

D 










- 


- 


C 






- 


0 








- 


0 






i 




- 


i 






- 


ak A 

1 1 1 1 








- 


1 1 1 1 


1 1 1 1 


1 1 1 1 


1 1 1 1 



FIG. 6: Euro and DAX futures. Factorial moments are displayed for 1, 2, 4, 10 or 20 bins. The global time window of 
16.5 hours {ribins = 1). This provides a time scale ranging from 1.6 hours for 20 bins till 16.5 hours for 1 bin. For the DAX 
future, the intermittent behavior is more pronounced than for the Euro. 

A final comment is in order. This analysis has been illustrated on the Euro futures. However, we have found that a 
similar intermittent behavior is observed on other futures, like the DAX futures (FDAX). For each series, the values 
of the like-sign and unlike-sign factorial moments of order 2 vary compared to results exhibited for the Euro future. 
However the rise of the F 2 as a function of the number of bins is an invariant property, with always and F 2 
equal to unity for a time window of order 16 hours. For some series, like FDAX, the intermittent behavior already 
prevails for a time scale of order 8 hours, with an increase of ^ 2 ^“'' from 1 {ribins = 1, At ~ 16 h) to 1.16 {ribins = 20, 
At ~ 1.6 h). This is shown in Fig. [0] This growth with ribins is faster than the Euro future. However, the key feature 
is always the same. For a sufficiently fine time scale, more precisely for At smaller than 4 to 8 hours depending of the 
financial product, intermittency is observed. 

Interestingly, we can use this formalism in order to derive some further statements on correlations of returns. From 
Eq. (HU, we get 

G{-l)=po 

which corresponds to the probability to have zero positive (resp. negative) returns in a given time window. This 
defines a gap in duration for positive (resp. negative) returns. Let us use Eq. (fTfl) to express G(— 1) in another way 
using cumulants Kq 

G{-1) = exp{-Ki + K2/2\ + ...) 

When we can neglect correlations within a large time window, we have shown K 2 = 0, then we conclude 

Po ~ exp(— ATi) = exp(— < >) = exp(— pAt) (24) 

where At represents the time scale (time window). This last expression is very simple and instructive. It states that 
the gap probability is exponentially suppressed with this time scale. This is illustrated in Fig. Qfor the finance case. 
We observe the distribution of events (probability distribution) as a function of the size of the gap. The gap is defined 
as the number of consecutive positive returns after a negative return, thus this is a gap in negative return. Inversely 
for a gap in positive return. This gap is given in number of time units for the financial time series considered. In Fig. 
[71 we display results for the Euro futures series over 10 years, sampled in two different time units. In both cases we 
observe effectively an exponential fall of the probability distribution as a function of the gap size. As illustrated also 
in Fig. [71 this exponential fall does not depend on the sampling. This confirms the prediction of Eq. (IMl) . 
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FIG. 7: Left: Euro futures sampled in 5min time units. Right: Euro futures sampled in 30min time units. The distribution of 
events (probability distribution) is displayed as a function of the size of the gap. The gap is defined as the number of consecutive 
positive returns after a negative return, thus this is a gap in negative return. Inversely for a gap in positive return. The Euro 
futures series over 10 years is sampled in two different time units (5min-Left and 30min-Right). In both cases, we observe 
effectively an exponential fall of the probability distribution as a function of the gap size. We have presented an exponential fit 
on top of each plot. On the Right plot (30 minutes sampling), we super-impose also the fit obtained on the 5 minutes sampling, 
in order to show that the exponential slope is identical in both cases. 



IV. FURTHER ANALYSIS OF CORRELATIONS 

With factorial moments, we have discussed how the deviation from Fq = 1 with the time scale At is associated to 
a broadening of the probability distribution of returns and then to dynamical fluctuations. 

Correlations between returns can be tested in other ways. One general methodology consists in estimating how 
a certain fluctuation measurement, labeled generically as F, scales with the size r of the time window considered. 
Specihc methods, such as the Hurst rescaled range analysis or the Detrended Fluctuation Analysis [13, differ 
basically on the choice of F{t). If the financial time series is uncorrelated one expects that F ~ as is the case 
for the standard Brownian motion. On the other hand, if F ^ with iL 7^ 1/2 one then says that the time series 
has long-term memory. The exponent FI is generally referred to as the Hurst exponent. 

To understand simply the interest of this exponent in describing correlations, we recall that the fractional Brownian 
motion (FBM) > 0| is a Gaussian process with zero mean and stationary increments whose variance and 

covariance can be written nj-H 



(25) 

< BH{s)BH{t) >t= i -I- -\t- , (26) 

where 0 < < 1 and < ■ >t denotes expected value. 

For H = 1/2 the process B}j{t) corresponds to the standard Brownian motion, in which case the increments 
X{t) = Bn/t -b 1) — Bn/t) are statistically independent and represent the usual white noise. For H ^ 1/2 the 
increments Xt display long-range correlation with 

< [X{t F T)X{t)] >t~ 2H{2H - 1)t 2^-2 for t ^ 00 . (27) 



where the average < . >t is made over t and the full time window of the analysis. Thus, if 1/2 < iL < 1 the 
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increments of the FBM are positively correlated and we say that the process Bnit) exhibits persistence. Likewise, 
for 0 < < 1/2 the increments are negatively correlated and the FBM is said to show anti-persistence 

Let us apply the discussion above on FBM Bnit) to real price series presented in previous sections. It is well known 
that the Hurst exponent can be extracted from financial price series using various techniques In doing so, we 

need to compute H from 



< x{t, rf >t=< [log[5'(t -k r)/S'(t)]]^ >t(x (28) 

in the range of validity of this expression. Obviously, we can write 

t'=t 



At this level, we can safely assume that the returns x{t') is identically distributed. However, the returns x{t') are 
not independent and thus the central limit theorem does not apply. In order to lead intuition, we shall anyhow make 
the approximation that we can sum up the sequence of x{t') in Eq. (12911 using the central limit theorem. Then, there 
are two cases (i) if the power tail of the distribution of returns follows v > 2, the variance is finite and x{t,r) 

converges to a Gaussian distribution; (ii) if z/ < 2, the sum in Eq. (|2^ tends to Levy stable distributions. As we 
have shown in previous sections for the futures under consideration in this letter, we stand in the case (i). Therefore, 
following the discussion above, we expect a Hurst exponent of order 1/2. 

In order to determine the Hurst exponent experimentally, we compute 



X ^ x/< \x{t + t) - x{t)\^ >t 

V< |a:(0P >i 



(30) 



Results are presented in Fig. [8]for the DAX and Euro futures. Then, the Hurst exponent is extracted with the relation 
F{t) oc . We get: H = 0.54 ± 0.04 for FDAX and H = 0.51 ± 0.03 for EC. These results are well compatible with 
the 1/2 value expected from the above discussion. They are also robust within the set of high liquid futures discussed 
in this letter. Let us note that the uncertainties on H are dominated by systematical effects. We have varied the 
time window for the calculation, the time interval and we have also divided the original data sets over 10 years on 
periods of 3 years. Also, we have modified the definition of the time average dehned in Eq. (13011 using an exponential 
smoothing in order to take into account the fact that the recent past is more important than the remote past. The 
variations in the determination of the Hurst exponents give the errors quoted above. 
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FIG. 8: Left: DAX futures. Right: Euro futures. In these plots F{t) = ^ Hurst exponent is extracted 

V<h(t)P>t 

with the relation F(r) oc r^. We get: F[ = 0.54 ± 0.04 for FDAX and H = 0.51 ± 0.03 for EC. See text for a discussion on the 
uncertainties. 



13 



Eq. (I30|) can be extended to higher moments of the distribution of returns as 



F,{t) 



yj< \x{t + r) 

V< >t 



(31) 



which leads to the definition of generalized Hurst exponent, namely -FIj(t) oc Then, processes with H(q) = H 

independent of q are called uni- fractal whereas processes where H{q) is not constant are referred to as multi- fractal 
When considering all systematic effects listed above, we do not get any experimental sensitivity to the q 
dependence in Eq. dSlD. 



V. CONCLUSION AND OUTLOOK 

We have presented experimental results based on a large set of time series on futures. When studying the probability 
distribution of log-returns, we have shown that a t-distribution with z/ ~ 3 gives a nice description of almost all data 
series for a time scale At below 1 hour. For At > 8 hours, the Gaussian regime is reached. A particular focus has been 
put on the DAX and Euro futures, which both verify these properties and thus exhibit a similar statistical texture. 
This appears to be a quite general result that stays robust on a large set of futures, but not on any data sets. In this 
sense, this is not universal. 

A technique using factorial moments defined on returns is described and similar results are obtained. Let us 
recall that factorial moments are convenient tools in nuclear physics to characterize the multiplicity distributions 
when phase-space resolution becomes small. In particular, correlations between particles lead to a broadening of the 
multiplicity distribution and to dynamical fluctuations. In this case, the factorial moments increase above 1 with 
decreasing resolution. This corresponds to what can be called intermittency. A similar property has been illustrated 
on financial price series. An intermittent behavior has been extracted for DAX and Euro futures, using moments of 
order 2 (F 2 ). This leads to perfectly consistent results with the standard distribution analysis. 

Let us note that from a fundamental point of view, there is no clear reason why DAX and Euro futures should 
present similar behavior with respect to their return distributions. Both are complex markets where many internal 
and external factors interact at each instant to determine the transaction price. These factors are certainly different 
for an index on a change parity (Euro) and an index on stocks (DAX). Thus, this is striking that we can identify 
universal statistical features in price fluctuations of these markets. This is really the advantage of micro-structure 
analysis to prompt unified approaches of different kinds of markets. 

Finally, we have discussed that the power law distribution of returns is well consistent with the information encoded 
into the Hurst exponent. We have obtained H = 0.54 ± 0.04 for DAX and H = 0.51 ± 0.03 for Euro futures. 

An immediate outlook concerns the use of this knowledge to unfold the sequence of returns into a more controlled 
sequence of trades. In Ref. [ 2 ^, we have proven that it is pos sible to design such unfolding procedure. Interestingly, 
the trade durations derived on the DAX and Euro futures correspond perfectly to the time scales derived in this 
letter. In fact, building a trading strategy is another consistency test of the knowledge we can get from the time series 
analysis. The broadening of the probability distribution of returns when decreasing the time scale and the associated 
dynamical fluctuations, as extracted from the factorial moments can serve as a guide in the unfolding procedure. 
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